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1. Introduction 

The search for the existence of the particle candidates comprising dark matter (DM) is still 
an open question. Many theories for DM involve extensions of the standard model (SM) such as 
supersymmetry (SUSY) [1] or theories with extra dimensions [2] with weakly interacting massive 
particles (WIMPs) that would annihilate or decay into standard model particles, most producing a 
continuum of y rays with energies up to the DM particle mass [3] or mono-energetic /-ray lines 
[4]. The indirect search for SM particles from astrophysical objects with a large inferred DM den¬ 
sity is a important complement to direct searches for DM interactions and accelerator production 
experiments. 

Attractive objects for indirect DM searches are nearby massive objects with high inferred DM 
density which are not known or expected to be /-ray emitters. The Galactic Center is expected 
to be the brightest source of / rays resulting from DM annihilation; however VHE emission ex¬ 
ists coincident with the position of SgrA* [5] and VHE emission along the Galactic Center ridge 
[6] [7], makes DM searches in the Galactic center a complicated, but not impossible, prospect. 
Other targets for indirect DM searches are Fermi-LKV targets without counterparts at other wave¬ 
lengths, which are potentially DM Galactic sub-halos [8][9] and galaxy clusters [10]. Einally, dwarf 
spheroidal galaxies (dSphs) are relatively close (~50 kpc), and have a low rate of active star for¬ 
mation, suggesting a low background from conventional astrophysical VHE processes [11]. This 
work will focus on the results of observations of dwarf galaxies with the Very Energetic Radiation 
Imaging Telescope Array System (VERITAS). 

2. VERITAS Data and Analysis 

VERITAS has devoted a significant portion (~15% in recent years of the total dark time ob¬ 
serving budget) to the observation of indirect DM targets, with most of that time on dSphs. This 
time is divided between the ‘classical’ dSphs, such as Draco and Ursa Minor (UMi), which have 
large stellar populations and therefore a well-constrained DM density profile, and the ‘ultra-faint’ 
dSphs, mostly Segue 1, which has a large inferred DM luminosity because of its close distance, but 
a lower stellar population and therefore a larger uncertainty in the DM luminosity. The inferred DM 
luminosity is referred to as the ‘J factor’ or ‘astrophysical factor’ and is defined as for annihilafion: 

Jann = J j Pil,^fdlda (2.1) 

where p(/,D) is fhe DM densify profile, D is fhe solid angle around fhe dSph and I is fhe disfance 
along fhe line of sighf. VERITAS divides ifs fime befween fhe classical and ulfra-fainf dSphs fo 
give a larger coverage in RA bands and as a safeguard in case any one parficular dSph falls ouf of 
favor of being a good DM fargef. As an example, fhe ulfra-fainf dSph Segue 2 was speculafed fo 
have a large J facfor because of ifs relafively close-by locafion (~35 kpc) [12] indicafing a good 
DM fargef, fhen additional sfudies led fo a re-evaluafion and a much lower J facfor. [13]. 

During fhe period from fhe beginning of four-felescope operafions in 2007 fo summer 2013, 
~230 hours of qualify dafa of dSphs have been faken. Qualify dafa for fhis work requires nominal 
hardware operafions of all four felescopes, clear' weafher, and low nighf-sky background (NSB) in 
order fo gef fhe lowesf possible energy fhreshold. The dafa span fhree VERITAS epochs: before 
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Dwarf 

Eive time 

logic-/ 

Significance 

r’95% 

-0-12 


[hrs] 

[GeV^cm^^] 

[rt] 

[lO^^^cm^^s^^] 

Segue 1 

92.0 

1 Q 4 - 1 - 0.3 

0.7 

0.34 

Ursa Minor 

60.4 

18.9+°| 

- 0.1 

0.37 

Draco 

49.8 

18.8±0.1 

- 1.0 

0.15 

Bodies 

14.0 

18.2±0.4 

- 1.0 

0.40 

Willman 1 

13.7 

N/A 

- 0.6 

0.39 


Table 1: Summary of dSph data used in this work. J factors are from [13] with an integration radius of 0.5° 
around each Dwarf. Dwarf galaxy detection significance is calculated from Li & Ma equation 17 [22] and 
the integral flux upper limits at 95% confidence flux is in units of 10^'^cm^^s^* above 300 GeV, assuming 
a spectral index of -2.4, using the Rolke method with bounded intervals [23]. 


and after the T1 relocation [14] and after the upgrade of the array with FPGA-based topological 
triggers [15] and camera pixels with higher quantum-efficiency photomultiplier tubes (PMTs) [16]. 
The five dSphs used in fhis sfudy are Segue 1, Ursa Minor (UMi), Draco, Bodies 1, and Willman 
1. The observafions and resulfs of fhe dafa analysis of fhese dSphs are summarized in Table 1. 

Dafa analysis mosfly uses fhe sfandard VERITAS techniques [17], wifh a few nofable ex¬ 
ceptions. In order fo consfrucf a viable background specfrum for each observation run (which is 
required by fhe even! weighfing mefhod described lafer in fhe fexf) fhe cosmic-ray background was 
subfracfed using an annulus cenfered around fhe felescope array fracking posifion as opposed fo fhe 
fargef locafion [18]. Also, in fhis work, looser cufs opfimized for sofl specfral sources were used, 
in order fo obfain fhe lowesf possible energy fhreshold. The combinafion of deep exposures and fhe 
lower energy fhreshold revealed background sysfemafics in fhe cameras, namely a gradienf along 
fhe direction of zenifh angle of observation and deficil regions in fhe camera due fo optically brighf 
sfars. The 3.8 magnifude sfar q-Leonis located 0.68° away from Segue 1 was of particular concern. 
The zenifh gradienf was corrected for by using a zenifh-dependenf accepfance function [19] [20]. 
The brighf sfars were corrected for using a 2D ellipfical fif fo each shower image, as opposed fo 
a momenf analysis [21]. This has fhe added benefif of inferpolafing over camera pixels fhaf are 
shuf off due fo high currenfs or have raised cleaning fhresholds due fo higher NSB levels [19], and 
improving fhe PSF, particularly af higher energies. 

3. Dark Matter Search 

To search for dark matter in all five dSphs simulfaneously, an evenf weighfing mefhod is em¬ 
ployed. The full defails of fhe mefhod are in [24]; if is only summarized here. If improves on 
sfandard lACT fechniques for DM searches by ufilizing fhe individual energies and spafial proper¬ 
ties of fhe individual evenfs, as well as fhe insfrumenfal and asfrophysical background properfies. 
Each even! in fhe region-of-inleresf (defined as fhe ON region in mosf lACT analysis) is assigned a 
weigh!, w,-, based on fhe dwarf field if originafes from, ifs reconsfrucfed energy E, ifs reconsfrucfed 
angular disfance away from fhe dwarf, 0 and dwarf EOV if originafed from, v. The fesf slafislic for 
fhe hypofhesis of / rays being produced by dark matter is fhe sum of fhe weighfs from all evenfs: 
T = ^w;. The index i is summed over all ON evenfs for all dSphs. The mosf opfimal form of fhe 
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weight function is 


w{v,E,e) =log 


biv,E,d) 


(3.1) 


where s{v,E,6) and b{v,E,6) are the total number of expected signal and background events 
with those properties, respectively. The number of expected background events is extracted by a 
spectrum generated by the OFF data for each data run using the Crescent method described in the 
previous section. The background is assumed to be isotropic within this region, so b{v,E,6) is 
proportional to sin(6)r/6. 

The expected signal s{v,E, 6) is determined by a convolution of the dark matter annihilation 
flux with the VERITAS instrument response 


s{v,E,e) 


dN{v,E,e) 

dEdO. 


dE2n ?,m{d)dd 


(3.2) 


where the number of events reconstructed with energy E and angular separation 0 is given by the 
convolution 


dN{v,E,e) 

dEdO. 


Ej JQ.1 


dEtdQ.t 


E{EtA) 

dEtdQ-t 


R{EAEtA) 


(3.3) 


where the subscript t denotes true energies and directions and the function R is the VERITAS 
response which is a function of the effective area, live time, instrument PSF, and energy dispersion 
(the probability E given Et). Effective area, PSE and energy dispersion are generated from y-ray 
simulations. The function E is the y-ray flux from dark matter annihilation: 


E{Et,dt) _ (av) dNy{E)dJ{d) 
dEtdQ-t dE dQ. 


where dNy{E)/dE is the spectrum of gamma-rays produced per annihilation, dJ{6)/dQ. is the DM 
density profile, M is the dark matter particle mass and (cjv) is the velocity-weighted annihilation 
cross section. The single annihilation spectrum used in this work is from PPP4 [25] and the DM 
density profile is from [13]. An example of the form of w, as a function of E and 0 for a dark 
matter particle of mass 10 TeV is shown in Eigure 1. 

In order to search for annihilation or set limits on the cross section, we need to compute a 
probability distribution for T, with the null hypothesis being that the observed test statistic, Tobs, 
is due to only background processes. The detection significance is defined as the probability that 
T is less than Tabs under the background-only hypothesis (equivalent to (av)=0). The description 
of computing the probability distributions is detailed in [24]. This distribution is computed for 
specified values of M and choices of DM annihilation channels. 

The search for annihilation is performed by measuring Tgbs and comparing this with the prob¬ 
ability distribution for T due to background. A search for a individual dwarf galaxy is performed 
by setting weights from all other dSphs to 0. The search for dark matter is performed as a function 
of mass and with the consideration of a 100% branching ratio to heavy quarks {bb), leptons (t+t^) 
and a two-photon final state. Eigure 2 shows the combined search for four of the five dSphs in this 
study, using the data summarized in Table 1. Willman 1 is not used because the DM profile could 
not be reliably determined [26]. Eigure 2 shows a non-detection of DM in the mass range from 100 
GeV to 100 TeV at a significance level not exceeding the 2a level above the background. 
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Figure 1: Example of the form of w, for a dark matter particle of mass of 10 TeV as a function of recon¬ 
structed energy E and reconstructed angular distance from the dwarf 0. 





Figure 2: Results of the search for dark matter annihilation of four of the dwarfs listed in Table 1 as a 
function of DM particle mass given annihilation into quarks {bb, left), leptons center) or directly into 

7 -rays ( 77 , right). Detection signihcance of the DM hypothesis dehned as the quantile of the probability 
distribution of the background-only model. Probability is converted to a sigma value by using the inverse 
CDF of a standard Gaussian. The dashed lines show the uncertainty in detection based on the DM density 
profile. The solid line is the median over the density profiles. 


4. Limits on DM Annihilation Cross-section 

With the absence of a DM signal, we compute the limits on the annihilation cross section. The 
95% confidence limits are generated by performing the hypothesis test at every value of the cross 
section. The (av)-space is divided in two regions where the hypothesis can and cannot be rejected 
at 95% confidence. The hypofhesis fesf is performed by asking, for a given value of (cJv), whefher 
fhe probabilify fhaf T < Toi,^ is less fhan 5%. The boundary befween fhe fwo regions represenfs fhe 
95% upper limif fo fhe cross secfion. 

Figure 3 shows fhe 95% confidence limif obfained wifh fhe 216 hours of dSph dafa. Each panel 
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Figure 3: Preliminary observed cross section exclusion limits as a function of DM particle mass with 95% 
confidence for annihilation channels bb (top left), (top right) and directly to 77 (bottom). The blue 

bands indicate the Iff systematic uncertainty in the dark matter profile. In the cases of bb and the 

cross section limits for the previously published 48 hour exposure VERITAS observations of Segue 1 are 
also shown [27]. The Gray band represents a range of generic values for the annihilation cross-section in the 
case of thermally produced dark matter. 


uses a 100% branching fraction into various Standard Model final states. The shaded band for the 
limits represents the 1 a systematic uncertainty that exists because of imperfect knowledge of the 
dSph density profiles. Figure 4 shows the median observed DM limit for the channel with 
la and 2a statistical uncertainties (the limit has a 68% and 95% chance to be within the red and 
green bands, respectively). Note that the observed limit does not exceed ±2a of the expectation 
limit for any value of the DM mass. 


5. Conclusions 

Presented here for the first time are the results of using the event weighting method for lACT 
analysis, using information from several objects to perform a combined search and compute cross 
section limits. No evidence of DM is found in 216 hours of combined dSph data, and limits of the 
annihilation cross section have been computed. The event weighting method has the added benefit 
of using the VERITAS PSF and individual event reconstruction position in addition to using the 
individual energies of each event, leading to more accurate cross section limits. The combination 
of softer cuts used for the analysis, roughly twice the exposure for Segue and the event weighting 
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Figure 4: Preliminary observed cross section exclusion limits as a function of DM particle mass with 95% 
confidence for DM annihilation to The red solid line indicated the ‘expected’ limit, generated from 

the background events. The red band and green bands represents the Iff and 2ff probability of the limit 
being within those regions, respectively. 

method leads to much more constraining cross section limits than our previous Segue 1 analysis 
[27] as shown in Figure 4. This is also the first time that a line search (DM annihilation directly to 
7 rays) has been performed with VERITAS data. 
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